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As devices to control spin currents using the spin-orbit interaction are proposed and implemented, 
it is important to understand the fluctuations that spin-orbit coupling can impose on transmission 
through a quantum dot. Using random matrix theory, we estimate the typical scale of transmitted 
charge and spin currents when a spin current is injected into a chaotic quantum dot with strong spin- 
orbit couphng. These results have implications for the functioning of the spin transistor proposed by 
Schhemann, Egues, and Loss. We use a density matrix formahsm appropriate for treating arbitrary 
input currents and indicate its connections to the widely used spin-conductance picture. We further 
consider the case of currents entangled between two leads, finding larger fluctuations. 
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I. INTRODUCTION 

There has been much recent progress in the creation 
and control of spin currents. There have been demon- 
strations and proposals for producing spin-polarized 
currents both witbiii^i^ii^i^ and without time reversal 
symmetry (TRS).— Recent progress in measuring and 
controlling the spin-orbit coupling in semiconductor 
heterostructure o^'^'"*^" promises to enable a range of spin- 
tronic applications relying on the spin-orbit interaction. 
As such devices are considered and developed, it is im- 
portant to understand the role of coherent mesoscopic 
fluctuations in these systems. In this paper, we consider 
the effects of injecting either a spin-polarized current or a 
pure spin current into a two-dimensional ballistic region 
with strong spin-orbit coupling and consider the scale of 
the fluctuations of charge and spin currents transmitted 
through such a device. 

For example, these effects could be important for 
the Schliemann-Egues-Loss spin field effect transistor 
(SFET) proposal." In such a SFET, spin-polarized elec- 
trons are injected into a region (e.g., a diffusive wire or a 
quantum dot) with spin-orbit coupling. In the "on" state 
of the device, the Rashbc^i^ and fc-linear Dresselhau o^^i^^ 
spin-orbit couplings are tuned to be equal, and the spin 
polarization does not decay as the electrons cross the 
region, but instead undergoes a controlled rotatiouiii 
In the "off" state, the Rashba and fc-linear Dresselhaus 
strengths are tuned to be different, and the spin polariza- 
tion is lost while traversing the region due to the random 
spin rotations experienced by electrons traversing differ- 
ent trajectories through the dot. Ideally, the on state has 
a fully spin-polarized current exiting the device and the 
off state has no spin polarization in the exit current. For 
coherent 2D quantum systems, however, the decay of the 
spin current in the off state relies on having a sufhciently 
large number of channels to average together. In the ID 
limit, with two ideal one- dimensional wires, each having 
only one propagating mode, a fully spin-polarized current 
injected into the first wire is a pure state, so the transmit- 
ted current must have a spin pointing in some direction; 
this fact implies that no reduction in spin-polarization is 



possible in the coherent ID limit. Other limitations to 
the SFET proposal have been simulated by Shafir et ali^ 
In this paper we discuss the general problem of co- 
herent propagation of currents through quantum dots, 
focusing on the relationship of incident to exit spin- 
polarization of the currents. For the case of 2D ballistic 
chaotic scattering regions with strong spin-orbit interac- 
tion, we use random matrix theory to give analytic re- 
sults for the expected values of spin-polarization in the 
exit currents. Once we can describe the ingoing current 
in terms of a density matrix, all of the conclusions will 
follow. Thus, the problem is generally broken into two 
parts: first, find the relevant input density matrix for 
the system of interest; second, propagate that density 
matrix to find the output currents and polarizations. We 
choose the density matrix formalism to describe the in- 
put currents to the quantum dots, as it is flexible enough 
to describe any current in the noninteracting system. As 
an important example, we describe how to construct the 
density matrices representing currents produced from po- 
tentials apphed to (possibly spin-spht) reservoirs. We 
go beyond this model and also consider injection of spin 
currents entangled between the two leads, finding larger 
fluctuations in this case. Similar work in a three-terminal 
geometry was considered in Ref. il6t . The case of unpo- 
larized input currents was considered in Ref. [^. 



II. SETUP 

We consider a quantum dot attached to two ideal leads 
through quantum point contacts (QPCs). There are A'', 
M open spin-degenerate channels in the left, right QPCs, 
respectively, and we let X = A^ -I- M. We take a basis 
for the propagating states in the ideal leads normalized 
to unit flux in each channel, as usual. We consider non- 
interacting spin 1/2 particles which are coherently scat- 
tered by the quantum dot, which we describe using an 
S-matrix. Given a density matrix w representing the 
current into the dot from the K channels, the output 
current is described by density matrix = SwS"^ . 

With K open channels, the S-matrix S can be repre- 
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sented by a 2K x 2K matrix of complex numbers. In 
systems with time reversal symmetry, it is convenient 
to consider S to he a, K x K matrix of 2 x 2 matrices. 
Any 2x2 matrix can be written as a linear combina- 
tion of the four Pauli matrices, but it is convenient to 
consider the basis {a'^ ,ia^ ,ia^ ,i(T^}, where the cr* are 
the Pauli spin matrices. In this basis, a 2 x 2 matrix 
q — q^a^ + iq ■ a, with g", (f G C, which is also called a 
quaternionJi Then q is defined to have a complex con- 
jugate q* = q°*tT° ■ a, dual = g°cT° -iq- a, and 
Hermitian conjugate q^ = q^* . The Hermitian conju- 
gate is the same as the standard Hermitian conjugate of 
a complex matrix, but the complex conjugate is not the 
same. For an S-matrix of quaternions, we define com- 
plex conjugate iS*)^J = iS^j)*, dual iS%j = (S^jO^, 
and Hermitian conjugate 5^ = S^*. This representation 
is convenient because for time reversal invariant systems, 
S = S^. The quaternion representation has the standard 
convention that tr(S') ~ J2i ^uj which is half of the trace 
of the equivalent complex matrix. 

A. Constructing w from chemical potentials 

Consider for the moment not two leads attached to the 
dot but K leads, each with one open channel and con- 
nected to its own reservoir with adiabatic, reflectionless 
contacts. Modeling the reservoirs as paramagnetic, each 
reservoir can be spin-split along its own quantization axis 
with each spin band separately in equilibrium, having its 
own chemical potential /ij^, where m e {1 . . . K} labels 
the channel and v G {0,x,y,z} indicates the charge and 
spin potentials. ^^'^^ There has been some confusion^° on 
the consistency of defining this chemical potential, so we 
give an example. If reservoir m is spin-split along axis x, 
then is the average chemical potential in the reser- 
voir, 2/x^ is the chemical potential difference between 
spin-up and spin-down electrons quantized along x, and 
f^nf — 0- general, if the quantization axis is h and 
the chemical potential difference along that axis is 2/1*^, 
then /i* — fj,^(h ■ i). Such spin-split chemical potentials 
can be realized, for example, by optical excitation in hct- 



erostructures, in an environment with inelastic relaxation 
much faster than spin relaxationji^iS^?^ 

We assume the leads have negligible spin-orbit cou- 
pling and spin relaxation, so there is a well-defined spin 
current in the leads. In the absence of inelastic processes, 
we can consider the current carried by particles with en- 
ergy e. For simplicity, we assume the number of open 
channels does not vary over the range of e considered 
here. Then the particle-currents flowing in from each 
channel are represented by the quaternion density ma- 
trix 

Wn,n{e) = <5„„,[/(6 - - a • flU'i^ - f^n)l (1) 

where /(e) is the Fermi function at temperature T, and 
we assume that /j,^ < max(T, A), where A is the mean 
orbital level spacing in the quantum dot without leads 
attached, and the prime indicates the derivative with re- 
spect to e. 

The charge current in the n^^ channel of particles with 
energy e is 

j°(6) = 2tr{P„[zi(6)-*™*(e)]}^, (2) 

where — e is the electron charge, h is Planck's constant, 
and P„ is the projection matrix onto the n'^ channel (i.e., 
{Pn)ab = SanSbn)- Similarly, the spin-current in the n^^ 
channel is 

j;(e) = 2tr{P„aM^i(e)-^™*(e)]}f . (3) 

We choose units in which e = ft. = 27r, so Eq. [3] can 
describe both charge and spin currents if we let a'^ be 
the identity. 

The currents are the physical objects in the system, 
and we note that the currents are unaffected by adding 
any multiple of the identity to u'(e), since w""* = SibS^ 
and S is unitary. We can thus use the density matrix to 
represent the currents, but we do not need to maintain 
tru" = 1 or even that w has positive eigenvalues. In 
the case where there are only two leads, we can subtract 
/(e - Ai") from w, giving 



w 



-f (6-/) <5/i"-/' (6-/i;)a.Arf 



I 

where = {^l + A*2)/2 and S^'^ = /ij* — /i2- Note that if Sfj,'^ < {T, A}), then we can represent the currents by 
dfi° = then the average chemical potential in both leads integrating over energy in the density matrix, giving 



is the same, so no net charge flows and w is traceless. 

If we consider an energy range in which the S-matrix w = | ^ + ■ Mi ^ ^\ 

does not vary (i.e., the linear response regime, where \ ' 1^2/ 



(5) 
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and 

j^,^2tiiP,-,a^w-w°^% (6) 

Spin-polarized injection from ferromagnetic contacts 
does not immediately map onto the chemical potential 
formalism. It is clear that if a ferromagnet is in equilib- 
rium with a wire, connected by adiabatic contacts, it will 
not produce a spin current in the wire, since adiabaticity 
requires that the lowest energy levels remain filled. For 
practical injection of spin-polarized currents from a ferro- 
magnet to a normal metal system, a tunnel barrier at the 
contact is the most common form of non-adiabaticityi^i^^ 

We can consider a situation where the ferromagnet in- 
jects into a semiconductor, which serves as the reservoir 
for a wire connected to our quantum dot. If we con- 
sider the case where the semiconductor has an energy 
relaxation time Tg much shorter than the spin relaxation 
time Tjs, then the spin-polarized current injected from the 
ferromagnet into the reservoir can relax to two indepen- 
dent distributions with a spin-split chemical potential. 
This is the same assumption used for optical excitation 
of spin-split chemical potentials. We can then use the 
formulation in terms of potentials as described above. 

The tunnel barrier at the ferromagnet introduces a 
second complication, as it implies that the ingoing cur- 
rent in the wire contains particles injected directly from 
the reservoir and also particles reflected from the scat- 
tering region and reflected back from the barrier. The 
input density matrix thus needs to be determined self- 
consistently, including the effects of both reflections. 
Such effects can be included systematically, by using the 
Poisson kernel^"* rather than the circular ensemble de- 
scribed below and also including the TRS-breaking ef- 
fects of the ferromagnetic scattering. For a sufficiently 
large reservoir in the semiconductor, this reflection can 
represent a small perturbation to the input currents, and 
the procedure described below will be a good approxima- 
tion. 



B. Connection to spin conductances 

We can write a generalized Biittiker-type conductance 
equation^^ 

jr = J2Gi^'t'i-^M,^r, (7) 

k,p 

where G'j^^ is the conductance from lead k to lead I and 
spin p to V and 2Mi is the number of modes, including 
spin, in lead /. The absence of equilibrium charge or 
spin currents (since there is no spin-orbit coupling in the 
leads) implies 

Y,Gr^-2Mi^Q. (8) 
fe 



Further, the conservation of charge current implies that 

^G?,- = 2A4J,o. (9) 
I 

Specializing to the case of two leads with N and M 
modes in the left and right leads with potentials fi'^, fj,'^, 
respectively, we can express G^^ simply in terms of the 
S-matrix. Setting /ij^ ~ and /i^ = S^a, Eq. [7] gives 

^•)^, = G^L- (10) 

Eq. El says that w = (^'^"^'^ ^^^"j ^ a"' Pl, and by Eq. [g] 

we have = 2tT{a'' PnSa^'PLS^) = C^^. Similarly, 
C'j^j^ = 2 tr {(j'' PrS a" PrS''). If the system is time re- 
versal invariant, then S = S^, which imposes some rela- 
tions between the different conductance matrix elements. 
Since tr{A-^) = tr(yl), we have the Onsager-like relations 

cr," = h'^hPGii (11) 

for k,l^ R,L where h"" = (1, -1, -1, -1). 

We thus see that we can express all of the G'^^' in terms 
of traces over appropriate density matrices multiplying S- 
matrices. We will consider the current in the right lead 
associated with the input density matrix w, defined as 

r^^2tr[a''Pn{SwS^ -w)l (12) 

which is proportional to the outgoing current in the right 
lead after injection represented by w, where the sign is 
chosen so that outgoing currents to the right are positive. 

C. Purity of w 

We will see that the coherence properties of the cur- 
rents are important, so it is interesting to consider when 
w represents a pure state. Ordinarily, density matri- 
ces are defined (with quaternion trace convention) so 
2 tr p = 1 , and p is pure if — p. In our open sys- 
tem, normalization is a choice, and we set 2trw — t, 
where t gives the total current incident on the dot. We 
can also add any multiple of the identity to w without 
affecting the physical currents. Taking both these factors 
into account, w represents a pure state only if there is a 
real number a such that 

/ w — al \^ w — al , , 

V2tr(w - al) ) ~ 2tr(w-al)' 

This condition implies that the K x K quaternion matrix 
w represents a pure state only if 

1. w'^ = tw, 

2. = 2K-l '^' 

3. w is invertible and 3a £ TZ such that 

1 _ to-[t-2a(K~l)]l 
^ ~ -a[t-a(2K-l)] ' 
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III. RANDOM MATRIX THEORY 

Though for any particular quantum dot it is difScuh 
to determine the fuU scattering matrix exactly, if there is 
a small number of open channels in the leads connected 
to the dot, mesoscopic fluctuations should produce an 
appreciable spin polarization in the exit current. We can 
understand this by considering that the current from one 
of the input channels has some probability to exit into 
each of the M exit channels after undergoing some spin 
rotation. In the chaotic strong spin-orbit limit, there is 
no correlation between the entry and transmitted spin 
polarizations. Though on average the transmitted spin 
polarization is zero, in any particular case there will still 
be some residual polarization in some direction in the exit 
lead. When there is only a small number of channels in 
the entrance and exit, these residual polarizations can be 
large. We will find the root mean square spin currents in 
the right lead by averaging over the ensemble of coherent 
cavities with strong spin-orbit coupling. These fluctua- 
tions are due to mesoscopic interference effects inside the 
quantum dots. We are primarily interested in the time- 
reversal invariant case, but we will present results valid 
with and without TRS. 

We consider coherent elastic scattering of noninteract- 
ing electrons with no spin-relaxation in the leads. We 
consider the chaotic limit for the quantum dot, in which 
the electron dwell time — 2'iTh/KA is much longer 
than the Thouless time rxh = Ld/vp, where Ld is a typi- 
cal linear distance across the dot, is the Fermi velocity, 
and A = 2TTh'^/mA is the mean orbital level spacing in 
the quantum dot, with m the effective mass and A the 
area of the dot. We further assume the strong spin-orbit 
limit, where the spin-orbit time Tso is much less than r^. 
We assume that all of the channels have perfect coupling 
into the quantum dot. 

We are interested in the properties of the current in the 
right lead. For an input density matrix w, in addition to 
j^, we define the outgoing current 



Jo 



2tr {aTRSwS'') 



(14) 



and the current due only to the input state 



so 



j^^ ^ 2 tr [a'' Prw) (15) 
Jout ~ Jin- The charge current is j° and the 



spin current is jw We define = jw ■ The polariza- 
tion of the current in the right lead is Pw = jw/jZ- ^ 
small number of parameters of the input current are suf- 
ficient to describe the effects of any ui in a two-terminal 
configuration. In particular, we define 



t = 2tru; 
C = 2tr(iy2) 

E" ^ 2tY{cr''PRwV) 



(16) 
(17) 
(18) 
(19) 

2tv{a-PRw"PRa-w''), (20) 

where superscript R is the quaternion dual, t is the total 
flux incident on the dot, C is a measure of the coherence 
of the current, gives the incident charge/spin current 
from the right lead, E'^ and are more measures of 
coherence. By adding a multiple of 1 to w , we can choose 
— 2tr{PRw) — 0, and all results below assume this 
choice. Note that if current is incident only from the left 
lead, then D" = E" = F" = 0. 

We take averages over the uniform ensemble of all 
S-matrices in the strong spin-orbit limit, either with 
TRS (called the circular symplectic ensemble - CSE) 
or without TRS (called the circular unitary ensemble - 
CUE).™s2i Such averaging is readily performed experi- 
mentally by small changes of the shape of a quantum 
dot^^; the root mean square (rms) fluctuations also give 
a typical value to be expected for any one chaotic dot. 
An external magnetic field can easily break TRS, mov- 
ing between these ensembles. A convenient formalism for 
performing such averages was worked out by Brouwer and 
Beenakker."^ From that work, we need two averages. In 
the quaternion representation, for /i = ti:{ASBS'^) for 
A, B constant K x K quaternion matrices. 



(/i) 



1 



- 2K-1 
1 

K 



(/l)cUE = -tl"^t^-^ 



2tr A ti B -ti (A'^B)] (21) 

(22) 



The other average we need is of /2 = 
tiiAS B S'f) ti- {AS BS'f) for A, B constant K x K 
quaternion matrices. We find^l 
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2/CSE 



2A 



^ ({K - l}|8[tr^]2[trB]2 + 2tr[A'^] tr[B^] + 4[tT{AB^)f - S,tr[A] tr[B] tr[AB^] - 2tT[AAB"B^]\ 
{^2[tTA]'^tT[B^]+tT[A'^][tTB]'^ -Atr[A] tT[A"B^]^AtT[B] tvlA^B"] 



4tY[A] ti[B] ti[AB'^] + tr[AB^AB^] + trlAAB'^B'^] 



}) 



(/2) 



CUE 



Ac/ 



[4i^ (tr A)2(tr B)^ + K tT{A^) tr{B'^) - tT{A^){tr Bf - (tr Af tr{B'^)\ , 

I 



(23) 
(24) 



where As = K{2K - 1){2K - 3) and Ku = KiAK"^ - 1). 

We consider the mean and fluctuations of j'^. Using 
Eq.Ell 



Jin 

(Jout) = 



2K~5s 



'2K -V 



(25) 
(26) 



where 6s — ^ for averages over the CSE and 5s = 
for averages over the CUE. The relevant fluctuations to 
study are of Aj^ = - (j^), which satisfy 



Jout 



(jout) 



(27) 



Using Eqs.[23]and[2l we find 



Jout 



= ^ {M5^n [UIt^{K -5s)- 2MC + A5sE° 



AW 



{K -5s){2MC + 25sD''^) (28) 



5s[E'>{2K-l)-F'^]], 



where A — A5, Ku for the CSE, CUE, respectively. We 
note that (j^^j^^^yg docs not depend on Z?, or F. 
Combining this result with Eq. [211 



Me{l + 5s) 



K - 5s/2 
2MC{K - 5s) - Mt^ 



2MC + A5sE^ 



(29) 



6s[{D'' 



, 2K^ -m- 

2K - 1 



-£;"(2A'-1)-F'']} 



Eq. [29] is the main result of this work, and we will now 
look at its implications in some special cases. First, an 
arbitrarily polarized current incident from the left lead, 
as can be readily created by optical methods. Second, 
a pure spin current uniformly distributed between the 
leads. Third, a pure state pure spin current, with entan- 
glement between the currents incident from each lead. 

Case 1: Spin-polarized current 

For any current incident exclusively from the left, the 
total current t and the parameter C are sufficient to de- 
scribe mean and rms currents in the right lead. We con- 
sider the input current represented by 

where s is the polarization magnitude and direction of the 
input spin current. Note that t can be positive, negative. 



or zero, depending on the direction of the charge current 
through the device. For |s| — the current is fully 
polarized. 

For the density matrix of Eq. [301 C = {t^ + s'^)/2N, 
and D = E = F = 0. Applying Eq. [211 the mean spin 
current in the right lead is zero and the average charge 
current is = 2tM/{2K - 5s). The reduct ion of 

(jZ) a-s TRS is broken {5s 0) is the signature of weak 
antilocalization i^^i^^i^^ The rms spin current in the right 
lead is 

_, M[{M-5s)t^ + iK-5s)s^] 
I ^A • ^ ^ 

The fluctuations in the charge current are 

(Ajf ) = f |[4M^ - 5sm - )^)\e + (1 - |).^| 

In the case of an unpolarized charge current (s = 0) 
with TRS, spin current in the exit lead is forbidden when 
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FIG. 1: (color online) For the fully spin-polarized current rep- 
resented by Eq. |30] with t = s = 1 and time reversal symme- 
try, comparison of numerical (symbols) and analytical (lines) 
results for the mean charge current (left), rms spin current 
(middle) and rms spin polarization (right) in the exit lead, 
where M (N) is the number of channels in the exit (entrance) 
lead. An average over 50 000 S-matrices from the CSE was 
performed for each data point. The lines are from Eqs. 1251 
1311 and |33l The right panel shows that the expected spin 
polarization in the exit lead is still appreciable, even for sev- 
eral open modes in each of the leads. Also shown are the 
equivalent CUE results with N = 2, showing that the rms 
spin polarization is nearly unchanged by breaking TRS in 
this case. 



M — 1 due to the combined effects of time reversal sym- 
metry and unitarity^ '^^i'^^ as can be seen in Eq. [3T] We 
can consider a pure spin current incident from the left by 
setting i = 0. In that case, we see that 



J'U 



M{N - (5s)s^ 
7VA 



(32) 



showing the scale of charge currents produced from the 
pure spin current. Similar effects have recently been 
proposed to measure the spin conductance in a three- 
terminal geometryji^ We note that (j°^)cse = if 
A'' = 1, showing that a pure spin-current incident from 
a single channel cannot produce a net charge current in 
the other channels. This is the time reversed statement 
of the theorem that with TRS a charge current cannot 
produce a spin-polarized current when A/ = 1. 

We can further consider the spin-polarization of the 
exit current, Pw = jw/ j^- It is clear that {pw) = 0, just 
as {jw) = 0, but there is some rms spin polarization of the 

exit current. If we approximate (p^) ~ (j^,^) / {iw) ' 
can use the above results to find 



{pl)^Z{K-5sl2f 



e{M-5s) + s\K~5s) 
M^MN 




2 4 6 
M 



8 2 4 6 8 
M 



FIG. 2: (color online) For the pure spin current represented 
by Eq. 1341 comparison of numerical (symbols) and analytical 
(lines) results for the rms charge current (left), mean spin 
current (middle) and rms spin current fluctuations (right) in 
the exit lead, where M (N) is the number of channels in the 
exit (entrance) lead. An average over 50 000 S-matrices from 
the CSE was performed for each data point. The lines are 
from Eqs. 1351 - 1381 The left panel shows that this pure spin 
current should still be expected to produce significant charge 
currents, with a nonmonotonic dependence on the number of 
open channels and M. Also shown are the equivalent CUE 
results with N = 2. 



To test this approximation, we found (p^) by numeri- 
cally averaging over the CSE. Matrices drawn from the 
CSE were chosen by diagonalizing matrices from the 
Gaussian Unitary Ensemble, as described in Ref. 1^. Re- 
sults are shown in Fig. [T] for the case t = s = 1, and 
it is clear that Eq. [33] agrees very well with the numer- 
ical results (right panel). The case shown in the figure 
is the relevant one for the Schliemann-Egues-Loss SFET, 
in which a fully polarized spin current is incident from 
one lead. In the off state, which relies on large spin-orbit 
coupling, the spin polarization in the exit lead is sup- 
posed to be zero. We see in Fig. [T] that even for several 
open channels in each lead, we expect to find an appre- 
ciable spin polarization in the output, limiting off-state 
function. 



Case 2: Pure spin current from both leads 

We consider a pure spin current incident from both 
leads, represented by the density matrix 



N 



2N 



2M 



(34) 



(33) 



This density matrix represents a spin current of +z in- 
cident from the left and a spin current of —z incident 
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from the right, which together are an incident pure spin 
current from left to right with polarization +z. In this 
case, t ^ 0, C = K/2MN, = (0,0,0,1), = 
(l/2Af, 0,0,0), and = (1, -1, -1, l)/2Af. Though 
the mean value of the charge current is zero, since it 
is as likely for the charge current to flow in as out, the 
spin current can produce a mean square charge current 



K 
1 



+ M 



21 



MNK 



(35) 



We note that when M = N = \, (jo)cse ~ showing 
that no charge current can be produced. This result is 
another implication of the theorem that, with time re- 
versal symmetry, a spin current incident in one channel 
cannot produce a charge current, combined with a sim- 
pler result that coherence and time reversal symmetry 
forbid spin-to-charge reflection in a single channel. 

The spin current in the right lead is a combination of 
the incident spin current, the reflected spin current from 
the right and the transmitted spin current from the left. 
Together, these give a mean spin current of 



(j-;) = (o,o,i-<5s 



1 



2K -r 



(36) 



Thus, with TRS, the spin current in the right lead is, on 
average, reduced from 1. In the case M = N = 1, this 
reduction removes 1/3 of the spin current that began in 
the lead. 

The fluctuations around the mean are 



K{K-6s){l-6sT!hi) 



NA 



K 



Ss 



{K -Ss/2)NA 
(M - 1)K^ 



M 



1 MI 
2M' ~ 



(37) 
(38) 

N 
K 



These results, along with confirming numerical simu- 
lations, are shown in Fig. [2l 



Case 3: Pure state pure spin current 

We consider entanglement between the currents in the 
two leads, which is beyond the standard chemical po- 
tential formulation of transport. In particular, consider 
a pure state spin current entangled between both leads, 
rather than the mixed state spin current of case 2. With 
M ^ N — 1. we consider 



1 



lay 



(39) 



This state has, as in case 2, a pure spin current +z in- 
cident from the left and a pure spin current —z incident 
from the right, but the off-diagonal terms of indicate 
that the currents are entangled. The density matrix for- 
malism easily allows consideration of such off-diagonal 



correlations between the channel currents. The entan- 
glement could be produced by passing a current through 
a beamsplitter produced from quantum dots^i^i^, feed- 
ing into the two channels or from spin injection by optical 
orientation using entangled photons. The density matrix 
represents a pure state by condition 3 of section III CI 
with a = —1/2. 

In this scenario, 1 = 0,0 = 3,0" = (0, 0, 0, 1), E" = 
(3/2,0,0,1), and F'' = (1, -1, -1, l)/2. This should be 
compared with case 2 in the AI = N = 1 limit, which is 
the same except C = 1 a.nd E" = (1/2,0,0,0). 

The most significant difference from case 2 is that co- 
herence between the channels allows a charge current 
to be produced, even when M = N = 1 with TRS. 
There is still no mean charge current, but the rms flue- 

tuations are (j°')c4e = 0-41, {fJ)cuE = 0-45. This 
rms charge current is much larger than the results of 
case 2, even away from M = N = 1, (see Fig. [5J left, 
and Eq. I35p indicating that the entangled spin current 
is better able to couple into charge current than is the 
incoherent spin current. We have normalized to have 
(jw) = (0:0,1 — Ss/3), as in case 2. We find fiuctua- 

■ 2 1/2 

tions around the mean of {Ajl^ )Jgg = (0.58, 0.58, 0.62), 

{AjlYJvE = (0.63,0.63,0.63). With TRS, the fiuctu- 
ations are larger along the polarization axis, but not 
markedly so. The total spin polarization fluctuations are 

(Aj-2^J,2= = 1-03, (Aj:')cue =.1-10 ^hich is larger 
than the mean current and equal in scale to the input 
current j?^, showing that coherence between the chan- 
nels signiflcantly enhances the mesoscopic fluctuations; 
this should be compared with Fig. [2] (right panel) . Such 
large fluctuations entail a signiflcant loss of knowledge of 
the quantization axis of the spin current, so the initially 
z-polarized current can exit polarized in many directions. 



IV. DISCUSSION 

Mesoscopic fluctuations of spin current on passing 
through a chaotic ballistic quantum dot can produce 
large fluctuations in spin-polarization, charge currents 
from pure spin currents, and spin currents from charge 
currents These predictions for mean and rms currents 
will be modified by dephasing and effects of the en- 
ergy dependence of the S-matrix. Dephasing processes 
can be readily added to this model using the third lead 
methodr^i^^i^ as detailed in Ref.[l. Dephasing generally 
reduces the fluctuations in charge and spin currents and 
also removes the symmetry that forbids charge or spin 
currents at certain values of M and N with TRS. 

If the ingoing current contains particles with energies 
varying over a large enough range, the energy depen- 
dence of the S-matrix must be considered as well. The 
S-matrix is generally correlated on the energy scale of the 
level broadening of the quantum dot eigenstates, approx- 
imately A' = AK/2 + 7^/2, where 'y^ is the dephasing 
ratci^i^^ If the incident particles have energies that differ 
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by a large amount compared to the level broadening A', 
as can happen at sufficiently large temperatures or 5)1'^ , 
then the mesoscopic fluctuations are suppressed, as there 
are effectively more open channels for particles passing 
through the dot. 

Mesoscopic fluctuations producing spin polarized 
exit currents could be important for operation of a 
Schliemann-Egues-Loss SFET. To avoid this impact on 
the off-state polarization, such a device should have many 
scattering regions in parallel or operate in a regime with 
sufficiently large temperature, bias, or dephasing so as to 
reduce these mesoscopic effects. 
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